5 



Sandpiles and absorbing-state phase 

transitions: 
recent results and open problems 

Miguel A. Murioz 1 , Ronald Dickman 2 , Romualdo Pastor-Satorras 
3 , Alessandro Vespignani 4 , and Stefano Zapperi 5 

1 Institute Carlos I for Theoretical and Computational Physics 
and Depto. de E. y Fisica de la Materia, Universidad de Granada, 18071 Granada, Spain. 
2 Departamento de Fisica, ICEx, Universidade Federal de Minas Gerais, Caixa Postal 702, 

30161-970 Belo Horizonte, MG, Brazil 
Dept. de Fisica i E. Nuclear, Universitat Politecnica de Catalunya, Campus Nord, Modul B4, 

08034 Barcelona, Spain 
4 The Abdus Salam Int. Centre for Theoretical Physics P.O. Box 586, 34100 Trieste, Italy 
INFM, Dipartimento di Fisica, E. Fermi, Universitd de Roma "La Sapienza" , P.le A. Moro 2, 

00185 Roma, Italy 



Abstract. We review some recent results on the relations between sandpiles and a 
class of absorbing state phase transitions. We use the concept of fixed energy sand- 
piles (FES), in which external driving and dissipation are absent. FES are shown to 
exhibit an absorbing state transition with critical properties coinciding with those of 
the corresponding sandpile model. We propose a set of Langevin equations capturing 
the relevant features of this transition. These equations characterize the universality 
class of systems with an infinite number of absorbing states and a static conserved 
field coupled to the order parameter. Different models in this class arc identified, and 
strong evidence is presented showing that the Manna sandpile, as well as some other 
stochastic sandpiles, belong in this universality class. Finally some open problems and 
questions are discussed. 



INTRODUCTION 



Slowly driven by the input of new ideas and concepts, avalanches of papers 
have in recent years seen the light in the field of self-organized criticality (SOC) 
[1-8]. This activity was certainly justified by the huge conceptual and practical 
implications of the self-organization mechanisms. From a theoretical point of view 



two main questions were soon posed: (i) what are the main ingredients able to 
drive "spontaneously" a system to its critical state without need of fine-tuning any 
control parameter?, and (ii) what are the links between self-organized critical states 
and standard critical ones appearing in equilibrium and non-equilibrium systems?. 
In order to answer the first of the two aforementioned issues, several paths to SOC 
have been identified (a recent review of them, as well as a summary on experiments 
on SOC can be found in [8]). Certainly the most celebrated ones among them are: 
the presence of two infinitely-separated time scales [6], and extremal dynamics [9], 
though these are certainly not the only ones [4,8]. 

For example, in sandpiles [1-3], the archetype of SOC systems, the two infinitely 
separated time scales correspond to a slow scale at which grains of sand are added 
to the system, and a fast one for the relaxation of unstable piles which leads to a 
rapid redistribution of sand, and eventually to boundary dissipation. It has been 
shown that it is only in the double limit of driving rate, h, and dissipation rate, 
e going to + with h/e — ► 0, that scale invariant behavior is generated up to the 
system-size cut-off [6,10]. Under the effect of the external drive, the system jumps 
among stable (absorbing) configurations via avalanche-like rearrangements. Some 
of the sandpiles we will refer to are: the original Bak-Tang-Wiesenfeld (BTW) 
sandpile [1], the Manna model [2], the Zhang model [3], the Oslo ricepile model 
[11], as well as some variations of them. 

Despite of some very remarkable rigorous results [4,12] and some renormalization 
group analyses [13], the question of whether the critical properties of sandpiles can 
be related to more standard critical phase transitions, or more specifically, whether 
the associated critical (avalanche) exponents [1,2,5] can be ascribed to any known 
universality class, remains open. 

A priori, these are not straightforward questions to answer; in particular, the 
presence of two time scales makes it difficult to apply standard techniques for the 
study of scale invariant systems to sandpiles. Another intrinsic difficulty lies in the 
fact, that exponents controlling the deviation from criticality cannot be measured 
in slowly driven sandpiles, posed by definition at their critical point. 

FIXED ENERGY SANDPILES AND ABSORBING 

STATES 

In order to circumvent these difficulties we have studied the the so-called fixed 
energy sandpiles (FES) [14-16]. These are cellular automata sharing the same 
microscopic dynamics as their corresponding standard sandpile counterparts but 
without external driving or dissipation. For each standard sandpile model we can 
define its FES version by keeping constant its total energy (number of sand grains). 
In this way, the initial energy density ( becomes a control parameter. FES are found 
to be critical only for a particular value C = Cc (which as we will show turns out to 
be identical to the stationary energy density of its driven-dissipative counterpart 



[16]). For energy densities below ( c the FES-sandpile falls into an absorbing state 
[17], while for values above that threshold activity is indefinitely sustained, i.e. 
the FES is in an active phase. Consequently, the phase transition in FES can be 
enviewed as an absorbing state phase transition [17,19] (with an infinite number 
of absorbing states). Observe that, switching on dissipation and driving, in the 
presence of activity ( > ( c and d(/dt < 0. In the absence of activity there is 
addition, but no loss of activity, so ( < ( c implies d(/dt > [18]. Evidently, the 
only possible stationary value for the density in the sandpile is ( c \ Not only ( c is 
common to a critical FES and its slowly-driven counterpart: also higher moments 
and correlation functions coincide in the large-size limit (as expected, given that the 
bulk dynamics is identical in both cases). We have performed extensive numerical 
analyses of the critical behavior of FES aimed at verifying the previous arguments 
and at determining a complete set of exponents. We refer the reader to [16] where 
the main ideas, simulations, analysis techniques, and results are presented (some 
of the obtained critical exponents are included in table 1). 

Let us remark that other attempts have been made to connect sandpiles with 
different types of phase transitions, as for instance: pinning-depinning transitions 
in disordered media [20], the voter model [21], and branching processes [22]. It is 
beyond the scope of this paper to review or discuss them. 

Generically systems exhibiting a transition into an absorbing state belong in the 
universality class of directed percolation (DP), the field theoretical (or Langevin 
equation) representation of which is the celebrated Reggeon field theory (RFT) 
[23,17]. This universality class has proven very robust with respect to the varia- 
tion of microscopic details, number of components, and number of absorbing states 
[24,25,17]. Only in the presence of some extra conservation law or additional sym- 
metry (as, for example, a Z 2 -symmetry between two equivalent absorbing states, 
which in some cases tantamount to parity conservation [26]), have the critical ex- 
ponents been found to differ from those of DP [26]. 

Following rather general principles of relevance and symmetry, a set of Langevin 
equations aimed at capturing the critical features of FES was proposed in [16]. Let 
us emphasize that this type of phenomenological approach has proven extremely 
useful, and provided sound results when applied to other systems with many ab- 
sorbing states (see [25]). This set of equations includes one for the evolution of a 
coarse-grained activity field, 0(x, t), and one for the energy (or background) field, 
-E(x, t) [16]. The activity field represents the local density of unstable sites, while 
the energy field stands for the total number of grains regardless of the local activity 
state. The proposed set of equations is: 

<%0(x, t) = fi(f)(x, t) - 60 2 (x, t) + L>V 2 0(x, t) + wE(x, t)0(x, t) + a y^x, t) 
d t E(x,t) = AV 2 0(x,t) (1) 

where /i,b, D,w,a, and A are constants, and r\ is a Gaussian white noise. In a 
nutshell, these equations express the fact that the transition for the activity is 
controlled by the same type of terms appearing in the RFT, i.e. the most rele- 



vant terms in standard absorbing-state phase transitions [17,23], plus an additional 
coupling between the activity field and a static conserved energy field. This extra 
term stems from the fact that creation of activity is locally fostered by the presence 
of a high density of the background field, and this background field (energy) is a 
conserved one. The extra conservation law is therefore the new relevant ingredi- 
ent with respect to RFT. Some other terms, consistent with the symmetries and 
conservation laws, could have been included in Eq.(l) but they all turn out to be 
irrelevant from a power counting analysis [16]. Let us also remark that the second 
equation in Eq.(l), being linear, can be integrated out, leading to a non-Markovian 
contribution and a quenched disordered linear term in the activity equation (a de- 
tailed analysis and discussion of these terms can be found in [16]). In what follows, 
we avoid performing such an integration to make explicit the presence of a con- 
servation law, and stress that this extra contribution is relevant at the RFT fixed 
point [16]. 

Which sandpiles is this theory intended to describe? Although in principle we 
would like to describe the original BTW sandpile this seems to be an unrealizable 
task, the reasons for this being manifold. The main reason, is that the BTW model 
being deterministic, has many hidden conservation laws and toppling invariants [4] 
which are almost impossible to implement in a coarse grained Langevin equation 
(for which a detailed book-keeping of symmetries and conservation laws is essential). 
Owing to the existence of these conservation laws, the FES version of the BTW 
exhibits strong non-ergodicity as shown and discussed in [16]; the stationary state 
depends strongly upon the initial configuration (and not only on the initial energy). 
From a different perspective, it has also been recently shown that BTW does not 
obey simple scaling [7]. 

Our approach is therefore intended to reproduce the scaling behavior of stochastic 
sandpiles with no extra conservation, symmetry, or spatial anisotropy [21], as for 
example: the Manna model [2], ricepiles [11], the stochastic version of the Zhang 
model [3,27], and other modifications of them. 

THE UNIVERSALITY CLASS 

Having proposed a set of Langevin equations for stochastic sandpiles, it would 
be an extremely desirable next step to derive specific predictions (i.e. values of 
the critical exponents and scaling functions), obtained after solving perturbatively 
the renormalization group associated with Eq.(l). Despite considerable efforts, 
unfortunately, so far we have not been able to complete satisfactorily such a chal- 
lenging task. Elucidating whether this is just due to technicalities, or has some 
more profound origin (concerning the theory renormalizability) still needs to be 
sorted out. 

Leaving aside for the moment field theoretical analysis of Eq.(l), we proceed with 
a slightly more indirect verification of the validity of our theory. Inspired by the 
previous set of equations Eq.(l) two of us have recently made the following con- 
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Manna 1.41(1) 3.36(1) 1.78(2) 1.76(2) 0.16(5) 0.78(2) 
DP 1.359(1) 3.507(1 ) 1.730(1) 1.901(5) 0.114(4) 0.730(4) 
TABLE 1. Critical exponents for spreading and steady state experiments 
in d = 2 and d = 3. Models: CTTP: Conserved threshold transfer pro- 
cess; CRD: Conserved reaction-diffusion model; Manna: Abelian Manna 
sandpile; DP: Directed percolation (for the sake of comparison). 



jecture: all stochastic models with an infinite number of absorbing states in which 
the order parameter evolution is coupled to a static conserved field define a unique 
universality class [28]. This conjecture has been supported by extensive simulations 
of the following different models fulfilling all the above- listed requirements [28,29]: 

• (i) a conserved threshold transfer process [30,17]; 

• (ii) a conserved lattice gas with repulsion of nearest neighbor particles [28]; 

• (iii) a reaction-diffusion model (CRD) with two species A and B [29] 

(see [28,29] for model definitions and details). All the measured critical exponents 
both for steady state and for analysis of spreading [17] coincide within error bars 
and are also the same as those measured in extensive numerical simulations of the 
FES Manna sandpile [16] (up-to-date exponent values in two and three dimensions 
are listed in table I; for further simulation details as well as model and exponents 
definitions see [16,28,29,17]). 

Furthermore, the slowly-driven versions of the previous models are self-organized 
[28,29], and the corresponding avalanche exponents are in good agreement with 
those of stochastic sandpiles [31,16]. Finally, all the scaling laws among steady- 
state, spreading and avalanche critical exponents [32] are satisfied within numerical 
accuracy (except for some well identified anomalies [16,28,29]). This shows rather 
unambiguously that the critical point of a slowly-driven sandpile coincides with 
that of its FES version, supports strongly the previously enunciated conjecture, 
and identifies stochastic sandpiles as belonging in the universality class of systems 
with many absorbing states and order parameter coupled to a static conserved field. 

Moreover, the conserved reaction-diffusion model studied in [29] has the noto- 
rious advantage of been susceptible to be mapped exactly using a a Fock-space 
representation and creation-annihilation operators [33] into an effective action, or 
equivalently into a set of Langevin equations [29] . Remarkably, up to naively irrel- 
evant terms, this set of equations coincides with Eq.(l). The additional irrelevant 
terms appearing in the exact mapping, might also be of some importance in lower 
dimensions (this issue is currently under study). 

CONCLUSIONS AND OPEN PROBLEMS 

Different stochastic sandpile models show critical properties that are fully com- 
patible with those of their corresponding FES versions. Using this equivalence 
stochastic sandpiles have been shown to belong in the universality class of systems 
exhibiting an absorbing-state phase transition with order parameter coupled to a 
static conserved density field. This result places sandpiles within the range of ap- 
plicability of standard tools for the study of scale invariance in phase transitions, 
and constitutes in our opinion a step forward in the rationalizing SOC. 

Finally let us enumerate some open questions and problems, on which we are 
presently working: 



• Direct numerical integration of the Langevin equations Eq.(l) using the 
method of [34]. This will provide an even more direct test for the validity 
of the Langevin equations. 

• Connection of the universality class discussed in this paper with that of the 
pinning-depinning transitions in disordered systems [20]. Numerical simula- 
tions show that these two classes are indeed the same; we are working on the 
theoretical connections between them [35]. 

• Renormalization of the field theory for FES, Eq. (1). 

• Clarification of some anomalies found in some time-dependent magnitudes in 
FES (like that affecting the exponent 9). 

• Extension of the numerical analysis to one-dimensional sandpiles. 
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